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Abstract 

The constraint equations for smooth [n + l]-dimensional (with n > 3) Rie- 
mannian or Lorentzian spaces satisfying the Einstein held equations are con¬ 
sidered. It is shown, regardless of the signature of the primary space, that the 
constraints can be put into the form of an evolutionary system comprised ei¬ 
ther by a first order symmetric hyperbolic system and a parabolic equation or, 
alternatively, by a symmetrizable hyperbolic system and a subsidiary algebraic 
relation. In both cases the (local) existence and uniqueness of solutions are also 
discussed. 


1 Introduction 

In solving the Cauchy problem in Einstein’s theory of gravity suitable initial data on a 
smooth n-dimensional manifold E, also refereed as an initial data surface, have to be 
chosen. The geometric part of this data can be represented by a pair of smooth fields 
( hij,Kij ) on E such that h %3 is a Riemannian metric while K VJ is a symmetric tensor 
held there. The constraint equations deduced from Einstein’s equation Gij — &ij = 0, 
with a source ^ki of vanishing divergence, are (for more details see |15] ) 

-e (n) R + (K j j) 2 - y-V - 2 e = 0 , (1.1) 

D j K j i -D i K i j -ep i = 0, (1.2) 

where e takes the value +1 for a Riemannian or —1 for a Lorentzian primary space, 

(n) 

whereas R and Di denote the scalar curvature and the covariant derivative operator 
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associated with h l3 , respectively. The source terms c and p, in (11.11) and (11.21) are 
given as 

c = n k n l &ki and p* = e h k in l @ki ■ (1.3) 

The above introduced universality might seem to be overwhelming especially as 
already in the Lorentzian case immediate conceptual issues emerge. These are rooted 
in the highly underdetermined character of the constraints. In particular, as it was 
emphasized in |2U] . the constraint equations impose only n+1 relations on the n(n+l) 
independent components of and K i3 , and the main issue is to single out from among 
these n(n + 1) independent components (or from their functions) those n 2 — 1 ones 
which can be freely specified, and those n + 1 ones which are determined by the 
constraints, and how to identify the ones which correspond to gauge transformations. 

A unique and almost exclusive approach applied as yet to solve the constraint 
equations is the conformal method The first essential step in devel¬ 

oping this method was made seven decades ago by Lichnerowicz mu- He constructed 
initial data for the n-body problem in general relativity by applying a conformally 
flat metric hij and vanishing extrinsic curvature Kjj. While applying the method 
proposed by Lichnerowicz one replaces the fields h %3 and the trace free part of K tJ by 

4 _ 

(p n ~ 2 hi 3 and <p” 2 K l3 , respectively, where ip is a smooth nowhere vanishing scalar field 
while hi 3 is an arbitrary Riemannian metric on £. Then, it was shown by Lichnerow¬ 
icz that a semilinear elliptic equation can be deduced from the Hamiltonian constraint 
for the conformal factor <p in the time symmetric case with vanishing extrinsic cur¬ 
vature. The second crucial step was made by York [211122] in the early seventies. He 
completed the method of Lichnerowicz by recasting the momentum constraint into 
a linear elliptic form, and putting thereby the full set of constraint equations into a 
coupled semilinear elliptic system. 

It is important to emphasize that by making use of the general method developed 
by Lichnerowicz and York one could, in principle, determine all the possible initial 
data configurations in general relativity. Nevertheless, as the momentum and Hamil¬ 
tonian constraints decouple if the trace of the extrinsic curvature K\ is constant—as 
well as, the momentum p, is York scaled (as described, e.g. in i) —most of the related 
developments are also based on this simplification. Therefore, almost all of the ex¬ 
istence, non-existence, or uniqueness theorems concerning the underlying semilinear 
elliptic system do also make use of these simplifications. For an excellent comprehen¬ 
sive recent review on the Cauchy problem and on the conformal method, in particular, 
on those results applicable in the non-constant mean curvature case, see Ref. [4]. 

The currently applied methods in determining initial data for binary black hole 
simulations are also based on the conformal approach and on the aforementioned 
simplifications. For instance, the most widely used method—based on the proposal 
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of Bowen and York [5]—assumes that the metric h %1 is conformally flat and that £ is 
maximal. As these assumptions are strictly stronger than the ones guaranteeing that 
the Hamiltonian and momentum constraints decouple we should not take by surprise 
that it is not possible to generate all the physically adequate initial data by making 
use of the corresponding method. Indeed, it was pointed out by Garat and Price in [6j 
(see also (9] for a more detailed argument), no conformally flat axisymmetric slicing of 
the Kerr black hole spacetime exists. Essentially this—along with the lack of a direct 
control of the physical parameters of the initial data specifications as a consequence 
of the implicit character of the conformal method, and along with the fact that one 
has to have both inner and outer boundary conditions and at the inner one we have 
much less support based on intuition than at the outer one—may be responsible for 
the non-negligible spurious gravitational wave content of the spacetimes evolved from 
Bowen-York type initial data specifications (for a discussion on this issue see, e.g., 
section 3.4.4. of p]). 

All the aforementioned observations underline that any alternative proposal al¬ 
lowing a convenient involvement of generic, i.e. non-conformally flat, metric h n] and 
generic type of extrinsic curvature K t] could be of obvious interest. In this paper 
these type of alternative methods will be proposed. In particular, it will be shown 
here that if the initial data surface £ is foliated by a smooth one-parameter family 
of homologous n — 1-dimensional surfaces then regardless whether the primary space 
is Riemannian or Lorentzian the 1 + n momentum constraint can be written as a 
first order symmetric hyperbolic system, whereas the Hamiltonian constraint may be 
viewed as a parabolic or an algebraic equation depending on our preferences on the 
separation of the controlled^ and freely specifiable variables in (11.11) and (11.21) . 

A key observation here is that as the constraint equations are highly underdeter- 
mined their character depend on the choice we make for the variables to be solved 
for. As a trivial example consider the equation 

dlf+d 2 y f+d x g-d 2 y g+h= 0 (1.4) 

which is elliptic, parabolic or algebraic depending on whether it is intended to be 
solved for the function /, g or h. respectively. 

Remarkably in both of the aforementioned cases the coupled Hamiltonian and 
momentum constraints can be solved as evolutionary systems. As it will be shown in 

*As it was shown in m the Hamiltonian constraint may also be put into the form of an elliptic 
equation for a conformal factor that can be defined by introducing the conformal structure on the 
foliating codimension-one surfaces (see also |T5] for more details on the conformal structure). Note, 
however, that the elliptic-hyperbolic systems yielded for the combination of this elliptic equation for 
the conformal factor and the hyperbolic system deduced from the momentum constraint is ill-posed 
M- Therefore, this elliptic formulation will not be considered in this paper. 
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section m if the Hamiltonian constraint is solved for the lapse of the foliation of £, and 
the level surfaces of the foliation are of positive mean curvature, then the Hamiltonian 
constraint can be seen to possess the form of a parabolic partial differential equation, 
and the coupled constraints form then a hyperbolic-parabolic system. Notably, the 
Hamiltonian constraint can also be solved algebraically for one of the components 
of Kjj. By eliminating this component, along with its derivatives tangent to the 
foliating surfaces, the momentum constraint, under suitable conditions, maintains its 
hyperbolic character. More precisely, it may be put into the form of a first order 
Friedrichs symmetrizable hyperbolic system. It is unique that in both cases a well- 
posed initial value problem can be seen to apply to the systems comprised either by a 
first order symmetric hyperbolic system and a parabolic equation or by a first order 
symmetrizable hyperbolic system and a subsidiary algebraic relation. 


2 The basic setup 

As it has already been indicated some restrictions concerning the topology of £ have to 
be imposed. In particular, we shall assume the existence of a smooth Morse function 
p : £ —» M that does not possess more than two isolated non-degenerate critical points 
on £.H This means that £ can be foliated by a one-parameter family of homologous 
codimension-one surfaces ZZ p —the latter denote the p = const level surfaces—such 
that (apart from possible critical points) the gradient Djp does not vanish on £. 

Note that the above assumptions—in spite of the involved technicalities—permit 
the involvement of high variety of topological spaces. For instance, all the prod¬ 
uct spaces Rxy are allowed, where the factor 5^ is a codimension-one manifold 
with arbitrary topology. This product structure guarantee that the “height function” 
determined by the factor M on £ will be a Morse function with no critical point. 

Notice that £ = M” belongs to the above set corresponding to the particular choice 
5? = M n_1 . Note, however, that M n may also be foliated—except at the origin—by 
codimension-one spheres § n_1 that corresponds to the Morse function p = 5T =1 (ah) 2 
with zero index at the origin in R n . Note also that the above assumptions permit the 
involvement of the n-dimensional sphere § n foliated by coclimension-one spheres S n_1 . 
This foliation corresponds to the choice of the height function p(x i,..., x n+ i) h->- x n+1 , 
where § n = {(aq,..., x n+1 ) G R n+1 | Y^=i( x i) 2 = ^ 2 } and the critical points are the 
north and south poles are represented by the points (0,..., 0, 1Z) and (0,..., 0, —7 Z) 
in M n+1 , respectively. 

In each of the above discussed cases the p = const level surfaces are, or they are 
supposed to be, orientable either compact and without boundary in £ or non-compact 

2 The critical points of a Morse function p (at which Dip = 0) are know to be isolated and non¬ 
degenerate in the sense that the Hessian of p is non-singular at those points. The index of a critical 
point is the number of the negative eigenvalues of the Hessian there m- 
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and infinite. H 

The rest of this section is to decompose—based on the use of the foliation 5^ p of 
£—the Riemannian metric h l3 and the primary extrinsic curvature . 

For instance, by making use of the unit norm vector field h l normal to the 5? p foli¬ 
ating surfaces, along with the metric 7 ^ induced on these level surfaces, the primary 
metric on £ can be decomposed as 

hij 7 ij T fliflj , ( 2 . 1 ) 

while the extrinsic curvature K, 3 of 5? p can be given as 

Kij = 7 Di fij = \ JZnlij • ( 2 - 2 ) 

Choose now a vector field p 1 satisfying the relation p l Dip = 1 throughout £. The 
unit norm vector field h l may also be decomposed as 

77 * = N~ 1 [p i — TV* ], (2.3) 

where the ‘lapse’ and ‘shift’ of p l are determined via the relations hi = N(dp)i and 
N l = 7 l j pi, with 7 l j = 8 l j — h l hj. 

In addition to the smooth function p : £ —y R. local coordinates (p,x 2 ,... ,x n ), 
adopted to the foliation 5? p and the vector field p\ can also be defined □ on subsets 
of £. This is done by choosing first (local) coordinates x A (with A — 2, 3, ..., n) on 
patches of one of the level surfaces =5^ of the foliation of £ and then by Lie dragging 
them along the integral curves of p l . 

The primary extrinsic curvature K tJ can also be decomposed as 

K.ij = k, hihj + [hi k j + hj kj] + K ;J , (2.4) 

where 

n = h k h l K kh kj = 7 fc j h! K kt and = j k a l j K u . (2.5) 

Notice that all the boldfaced quantities may be considered as if they lived on the 
p = const surfaces. We shall refer to the traces 

K\ = f l K kl and K.\ = f l K kl , (2.6) 

3 The 5f p level surfaces could also be allowed to be compact with boundary in S. Note that in this 
case—which, for conveniences, is left out of the scopes of the present paper—both the hyperbolic- 
parabolic and the symmetrizable hyperbolic systems can be solved (locally) as an initial boundary 
value problem such that the boundary values have to be fixed at the boundaries of the S^ p level 
surfaces in E. (For useful hints on the solvability of the parabolic problem in the pertinent case see, 
e.g. pages 25-42 in [7].) 

4 In these local coordinates the spatial indices of the pull backs of tensorial expressions to the 
level surfaces will be indicated by uppercase Latin indices and they always take the values 2,3,..., n. 
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and also to the trace free part of defined as 

K« = K„- Jt7«kV (2.7) 

Note finally that the n (n + 1) independent components of the original pair of 
variables (. hij,Kij ) may then be equivalently represented by the components of the 
variables in the septuple (N, N\ 7 ^; k, k,;, K |, Kj j). 


3 The constraints in new dress 


By making use of the above decompositions, along with a sensible adaptation of 
(A.9)-(A.10) and (A.12)-(A.13) of [13], the momentum constraint can be seen to be 
equivalent to the system 

^nk, - Di{K l i) ~ DiK + D l Ku + (K\) k, + K„ - e p, f, = 0 (3.1) 

^n(K \) - D% - K (k\) + K kl K kl + 2 h l ki + eP) h l = 0 , (3.2) 


where Di denotes the covariant derivative operator associated with 7 ^, and hk = 
n l Dihk = —Dk( In IV). 

Notice that only the source terms—the last term on the l.h.s. of (13.11) and (13.21) . 
respectively—pick up an e factor which means that the principal parts of the above 
equations are not affected by the signature of the metric of the primary space. 

Remarkably, 73 ^ A” times of (13.11) and N times of (13.21) when writing them out 
in (local) coordinates (p, x 2 ,..., x n ), adopted to the foliation S? p and the vector held 
p 1 , can be seen to take the form 


0 


d 0 + 


f-^N K *f AB -N j AK \ 

\ -N*i BK —N k J 




0. (3.3) 


Notice that the coefficients of d p and dx are symmetric, and that the coefficient of d p 
is also positive definite. Thereby (13.31) possesses the form of a first order symmetric 
hyperbolic system 

A {p) d p u + A {k) d K u + B = 0 (3.4) 

(with K = 2 ,,n) with the vector valued variable 

u = (k B .K E E ) T , (3.5) 


and such that the ‘radial coordinate’ p plays the role of ‘time’. 
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Notice also that the characteristic cone of (13.31) — apart from the S? p codimension- 
one surfaces (h*£j = 0 on them) foliating £—reads as 

[f-(n-l)nW]^- = 0. (3.6) 

It is also remarkable that all the coefficients and the source terms in (13.31) may be 
considered as if they were defined on the level surfaces 5? p exclusively. 

In turning to the Hamiltonian constraint recall first that in virtue of (A.l) of (13) 
for the scalar curvature^ of h, : j 

(n) R = R- {2 Jf A (K l t) + ( K\f + K kl K kl + 2 N' 1 D z Aiv} (3.7) 

holds, where R denotes the scalar curvature associated with 7 y. Combining this 
with the decomposition (12.41) of the primary extrinsic curvature K tJ the Hamiltonian 
constraint (II. 1 |) can be put into the form 

-e R + e {2 ^n(K\) + {K\) 2 + K kl K kl + 2 N 1 D l DiN^ 

+ 2 « {K\) + (K ^) 2 - 2 k% - K kl K kl - 2 e = 0 . (3.8) 

4 The evolutionary systems 

The system comprised by (13.ID - (13.2D and (13.81) has to be solved simultaneously. 
This section is to discuss two cases—each corresponding to a specific choice for 
the controlled variables in (j3. 8 |) — such that both of these lead to an evolutionary 
setup. In particular, if (I3.8j) is solved for the lapse N the corresponding equation 
becomes parabolic, and the coupled equations fl3.1l) - fl3.2D and (13.8p determine then 
a hyperbolic-parabolic system to which an initial value problem can be setup. Al¬ 
ternatively, (13.81) can also be solved algebraically for k, and the corresponding set 
of equations (13.ID - (13.21) and (13. 8 D turn to be a hyperbolic-algebraic system. Notably, 
by eliminating k and its derivatives from (13.ID - (13.21) . under suitable conditions, a 
Friedrichs symmetrizable hyperbolic system can be deduced from them to which a 
well-posed initial value problem known to apply. 

The following subsections are to deduce these hyperbolic-parabolic and symmetriz¬ 
able hyperbolic systems and to identify those conditions which guarantee the existence 
and uniqueness of local solutions to the pertinent equations. 

5 Here, as the metric hij on S was assumed to be Riemannian, regardless whether the metric of 
the primary space was of Euclidean or Lorentzian signature, the value of e—that applies in a sensible 
adaptation of (A.l) of (13] —is +1. 
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4.1 The hyperbolic-parabolic system 

This subsection is to show that if (13.81) is solved for the lapse IV, and if the S? p level 
surfaces of the foliation of E are of positive mean curvature, a parabolic problem 
manifest itself, and that the coupled equations (13.ID - (13.21) and (13.8D give rise to a 
locally well-posed hyperbolic-parabolic system. 

To see that this is indeed the case recall first that 

k„ = i % = Af - 1 !- d (i y, ] = i v-'ftq , (4,i) 

and that 

k\ = kij = N~ l [\ f j 2? p % - DjN j } = N-'k , (4.2) 

where K %] and K = 7 3 K,j were introduced only for notational conveniences. It is 
also straightforward to see that 

^n{k\) = -N~ 3 k [ (d p N) - (N'DtN)} + N~ 2 [{d p k) - (IV* A k)} . (4.3) 

By substituting the above relations into (13.8D we get 

—e R + e {—2 N~ 3 K[ (d p N) - (N l DiN) ] + 2 N~ 2 [ (d p K) - (IV* A A') ] 

+ N~ 2 K + N~ 2 K kl K kl + 2IV' 1 A Aiv} 
+ 2 « (K*,) + Af (K \f - 2 k*k) - K kl K kl - 2 c = 0 , (4.4) 
which by applying the notation 

4 = 2 [ (dpk) - IV* (A k) ] + K + K kl K kl (4.5) 

B= — R + e [2 n (K*/) + Af ( k M 2 “ 2 k*k) — K w K fc * « 2 e ] (4.6) 

can be put into the form 

2K[(d p N) - N l {DiN)] = 2N 2 (D l D l N) + AN + BN 3 , (4.7) 

~k 

which is a Bernoulli-type parabolic partial differential equation provided that K , 
or equivalently—as, it will be shown below, N can be guaranteed to be positive— 
K\ = N~ l K is positive throughout E, i.e. 5? p is a positive mean curvature foliation 
of E. Accordingly, if this happens equation (14.71) is uniformly parabolic such that p 
plays the role of “time”. 

Note that, in virtue of (14.21) . 

k = \Y>^ lj -D,N\ (4.8) 
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and that both terms on the right hand side can be evaluated once N l and 7 ^ are 
known on £. Since these variables will be seen to belong to the freely specifiable 
part of the data, as listed in (14.91) . they can be chosen so that K is guaranteed to 
be positive throughout £. This verifies then that equation (14.71) is indeed uniformly 
parabolic. If, in addition, smooth positive initial data 0 N is given on a level surface 
p = p 0 of the foliation of S, then standard results guarantee the existence of a unique 
positive solution for N, with N\^ = 0 N, on a short interval [p 0 ,Po + e) for some 

e > 0 , i.e. in a one-sided neighborhood S^[ P0iP0+e ) of S^ Po in £. 

Note that as the fundamental characters of the principal parts of (13.ID - (13.21) and 
(13.81) . respectively, is not affected by the coupling of these first order symmetric hyper¬ 
bolic and parabolic systems the following short time existence and uniqueness results 
apply to this coupled system: 

Theorem 4.1 Suppose that a smooth choice had been made for the freely specifiable 
variables 

N\^ i:j ;«,Ky ; c,p* (4.9) 

★ A . 

so that K > 0 throughout £. Assume that smooth initial data ( 0 iV, 0 kj, oK 1 1 ), with 
0 N > 0, had also been chosen to the hyperbolic-parabolic system on one of the level 
surfaces SA po in £. Then, on a short interval [po, Po + £)> f or some £ > 0, there exists 
a unique smooth solution (TV, kj, K Z /) ; with N > 0, to the hyperbolic-parabolic system 
comprised by (EJ]H22j) and UTM such that N\^ pq = 0 iV, k* | ^ = 0 kj,K z z |^ o = 
oKV The fields htj and Kij which can be built up from this solution and the freely 
specifiable part of the data satisfy the full constraint system \l.l\) - \1.2i) in a one-sided 
neighborhood y , [ P0tP0+£ ' ) of y po in'E. 

Clearly, it would be of obvious interest to know whether by imposing certain 
conditions on the coefficients and on the initial data the existence of unique solutions 
to this coupled hyperbolic-parabolic on the global time interval [po, 00 ) could also 
be shown. It would be even more important to identify those conditions which could 
guarantee that the initial data ( h tJ , K^) —that could be recovered from these solutions 
and from the freely specifiable part of the data—was of asymptotically Euclidean. The 
study of all these issues is out of the scope of the present paper and they are left for 
future investigations. However, as for certain choice of coefficients in (14.71) and the 
initial data 0 N on S^ po classical solutions are known to blow up in finite “time” it is 
of obvious interest to know those conditions which guarantee the existence of global 
in the time, [po, 00 ), solutions to (14.71) . 

In identifying these conditions we shall apply an argument adapted from Bartnik 
[ 2 ] (see also [T 8 ] CEH]). Notably many details of Bartnik’s proofs developed to study 
the existence of quasi-spherical foliations in asymptotically flat spacetimes do also 
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apply in the generic setup used in this paper. For instance, consider a function / on 
£ and, following Bartnik’s notation in [2], define f* and as 

f(p) =sup{f(p,x 2 ,...,x n )}, f,(p) = mf{f{p,x 2 ,...,x n )}. (4.10) 

s* P 

Substituting then N in (14.711 by applying the relation w = N~ 2 we get 

dpW + N^N 1 ^^) = ~2K~ 1 N~ 1 (D l D l N)-aw-b, (4.11) 

★ ★ 

where a = A/K and b = B/K. Notice that in the “gradient” and “Laplacian” terms 
the substitution was not performed as the parabolic maximum principle, along with 
the fact that N has a maximum where w has a minimum and, similarly, N has a 
minimum where w has a maximum , implies that at the maximum of N the relation 


d p w * > —a* w* — b* 

holds. This immediately gives us a lower bound 


u>* > IF* = exp 

l 

1 

2 

a 

* 

cn 

j 

W *\^ P0 ~ [ & * ex P 

fP' 

/ a*dp" 

dp'] 


L J po J 


k '-'PO 

j PO 

J 


By a completely analogous argument we also get an upper bound 


w * <W* = exp 

— / a* dp' 

\ 

W *\^P0 - [ b * eX V 

fP' 

/ a*dp" 

dp'] 


- J PO 


K J pa 

J pa 

J 


(4.12) 


(4.13) 


(4.14) 


Based on (14.1411 . and following the argument in |18| . we may conclude, on the one 
hand side, that w is bounded from above, w < IF*, which implies that N is positive 

_i _i 

as then 0 < ( W *) 2 < w 2 = iV. If, in addition, 


tc 


sup 

p£[po,°o) 


f fP 

fp' 

1 

/ b* exp 
( Ppo 

/ a* dp" 
j PO 

dp' 


(4.15) 


is positive and finite, i.e. 0 < /C < oo and if for the initial data w^\^ P0 the inequality 
> /C holds then it follows from (14.131) that w is also bounded from below by 

„ _i _i 

IF* > 0. This implies, in turn, that N has to be bounded, i.e. N = w 2 < IF* 2 < 00 
provided that for the initial data 0 N the inequality 0 N~ 2 = W \^ P0 > w * |. 5 % > /C, or 
equivalently qN < l/VK,, holds. Remarkably, if B , given by (14.61) . is negative then 
/C < 0 from which it follows that for any positive choice of the initial data qN to (14.71) 
IF* is positive and, in turn, the relation N < 00 follows. 

By combining all the above observations we get 
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Theorem 4.2 Suppose that all the coefficients in U. 71) are smooth and that the freely 

~k 

specifiable part of the data was chosen such that K is positive throughout E. Choose 
o N to be a smooth positive function on the level surface p = Po such that 0 N < 1/ VlC 
if 0 < /C < oo, or to be arbitrary if 1C <0. Then flA 7| ) has a unique smooth global in 
the time, [po, oo), classical solution such that 0 < N < oo, and that N\& p = qN. 


We close this subsection by comparing the assumptions used in our setup and 
those applied by Bartnik’s in [2] and by Smith and Weinstein in [TS1 (T9]. To do so 
recall first that in all of these works, as opposed to the setup proposed in this paper, 
attention was restricted exclusively to the time symmetric case with vanishing K l0 . 
Furthermore, in ®mm the foliating T/ p surfaces were assumed to be topological 
two-spheres whereas in our setup they were allowed to be generic as far as their mean 
curvature is positive. Moreover, the metric 7 ^ in Bartnik’s quasi-spherical foliations 
was supposed to possess the form 7 ^ = r 2 7 ^, where °f tJ stands for the unit two sphere 
metric and the function r : E —> M determines the level surfaces of the foliation of E. 
This setup was slightly generalized by Smith and Weinstein in [ITS| ITU] by assuming 
that 7 ij to possess the form 7 ^ = e 2v r 2 jij with some smooth function v : E —> M 
though they also required that the r = conts level surfaces have area 47 rr 2 . Both of 
these foliations are quasi-convex in the sense that the Gauss and mean curvatures of 
the level surfaces of the foliation SC p are positive, and the quasi-spherical foliation of 
Bartnik can be recovered by setting v — 0 in [T 8 l T9] . 


4.2 The symmetrizable hyperbolic system 

As mentioned above (13.811 may also be solved for n. Since (13.IK involves the “angular” 
derivatives of k as a first step in demonstrating the viability of the proposed method 
is to eliminate these derivatives and then to show the solvability of the yielded system. 

Once the derivatives Din are eliminated from (13.111 the reduced system can be 
seen to be comprised by (13.2H and 

+ (K^-^k A(K z 0 - 2k'A:k/] + (2K z i )- 1 A«o (4.16) 

+{K\) k f T [ k - At (K'j) ]hi-h l K u + D l K u - e pn l i = 0, 

where k, and kq are supposed to be replaced by the algebraic expressions 

^ = (2K i z )- 1 [2k'k i - g (K^) 2 - K 0 ] (4-17) 

f 7? 1 O O I I 

Ko = -f R - K kl K kl -2c. (4.18) 

Note that (j4.17(1 . with the substitution of (14.1811 . represents the algebraic solution of 

A. /s . O 

(13. 8 H for k, and that kq depends only on the quantities N,N l ,^ij and K^ which, in 
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the present case, will be seen to comprise the freely specifiable part of the initial on 
E. 

Note that the above outlined process guarantee that solutions n, kj, K \ to the 
system comprised by (14. 1611 and (13.2lh along with the freely specifiable part of the 

/V /s . O _ _ 

data IV, N l , 7 ^-. K^-. in virtue of the relations (12. IK and (12.41) determine solutions to 
the full constraint system (ll.ll) - (ll. 2 l) . 

Notably, when writing out (14.161) and (13.21) in local coordinates (p, x 2 ,..., x n ), 
adopted to the foliation 5A p and the vector field p®, they take the form 

dp\\ + B^ d K u + C = 0 (4.19) 

for the vector valued variable u = (kg, K £ g) T , where the coefficients B (R •' of 8k, with 
K = 2, 3,..., n, can be given as 


B\f = —N n 5ij - N^SinSj^ + 


rK 


KL 


kN 


Si(K-l) Sjn 2 


Nk 


L P LM 


K E , 


£ Lilvl r r 

7 Ci(A'-l) 0(j+l)M 5 
(4.20) 

and where the lowercase and uppercase Latin indices take the value 1,2,... ,n and 
2,3,..., n, respectively^ 


A system of the form of (I4.19p is known to be Friedrichs symmetrizable (see, 
e.g. j3j E, EZ]) if a common symmetrizer of the coefficients £? (A ) exists. A bilinear 
form, T-Lij, is called to be a common symmetrizer to (14.191) if it is positive definite and 
is such that for any value of K — 2,..., n the product TLij B^j is symmetric. It may 
be verified by a direct calculation that the matrix of the form 


Tiij 


Sin Sjn 


2 

- kp [ Si n S(j+1)Q + %+l)Q Sj n ] (4-21) 

tx> 

d-2 [ ^ (kp 1 PQ %+i)q) (ks 7 5T %+i)t) - « • K E e 7 PQ %+i)p %+i)q ] 


is indeed a common symmetrizer to (14.191) in a subset of E where neither n nor 
vanishes and in which they are of opposite sign. This observation is summarized in 
the following lemma. 


Lemma 4.1 Assume that n ■ K \ < 0 holds in a subset of E. Then equations \f.lb\ ) 
and \3. 6 A) comprise a first order symmetrizable hyperbolic system there. 

6 As the lowercase and uppercase Latin indices take the value 1,2,... ,n and 2,3,... ,n, respec¬ 
tively, in evaluating a term of the form here and also in (14. 2111 . one needs to refer to the 

(n + l)-dimensional Kronecker delta. 
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It is also informative to inspect the characteristic directions = (£ p ,£ 2 , • • • , £n) 
which are the roots of the essential part of the characteristic polynomial 

K\ (n%) 2 - 2 (n%) (i KL k K £ L ) + k ( 7 * L ZkZl) = 0 (4.22) 

of (14.191) . In virtue of (12.31) and (14. 22[) . for the p-component of a characteristic direc¬ 
tion the relation 


= n k £ k + 


N 


(*/ KL k K £ L ) ± \J(*f KL k K £ L ) 2 - k • K e e ( 7 kl £ k £l) 


(4.23) 


holds. Notice that in (14.23f) the inequality k ■ K \ < 0 comes into play again. 

The above discussion indicates that it is of crucial importance to know if this sign 
condition may be guaranteed to hold. 

To get some hints about the viability of the proposed setup note first that, in 
virtue of (12.11) and (12.41) . the relation 

K e e = k + K\ (4.24) 


holds. It follows then that for any maximal ( K e e = 0) but non-time symmetric 
(. K e f 0 ) initial data surface the signs of n and have to be opposite except 
where they both vanish simultaneously. 

As it would be preferable to involve more generic initial data surfaces than the 
maximal ones it is important that by choosing the freely specifiable part of data on 
E, along with the initial data okj and oK^ for the system (14.161) and (13. 2 j) . properly 
the inequality k ■ K \ < 0 can always be guaranteed to hold locally. To see this note 
first that the initial data ok; and 0 K^ for (14.161) and (13. 2 p is freely specifiable at one of 
the level surfaces of the foliation denoted by 5^ Po . Using this freedom choose oK^ to 
be strictly non-zero there. This, implies then that for any choice of a continuous held 
K. l i on E such that K. l i\^ po = qK^ the held does not vanish in a sufficiently small 
neighborhood of S^ po either. Similarly, by choosing the initial data for 0 k,;, along with 
the freely specihable data on E, properly the term k 0 — 2 k^ can also be guaranteed 
to be strictly positive on 5 ^’ po , and, in turn, in a sufficiently small neighborhood of 
y po in E. 

Viewing then (13.81) as quadratic polynomial 

(K^) 2 + k (K' t ) + [ - 2 k'k, ] = 0 (4.25) 

in K l i, and applying Viete’s formula relevant for the product of the roots (K^)! and 
(K l i) 2 we get 

(K i ;) 1 (K i ;) 2 = Ko - 2k'k;>0. (4.26) 
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This verifies then that the two roots are real and they are both either negative or 
positive in a neighborhood of S^ po . In virtue of the other Viete’s formula, pertinent 
for the sum of the roots, we have that 

(K^) 1 + (K z ; ) 2 = (4-27) 

which implies then that k and (K^) are of opposite sign at least in a sufficiently small 
neighborhood of S^ po in E. 

As for symmetrizable hyperbolic systems the initial value problem is known to be 
well-posed in the smooth setting, in virtue of the above observations, we immediately 
have the following: 


Theorem 4.3 Assume that the freely specifiable part of data 

N, N\ jij ; Kij ; e, pi (4.28) 

on E, and the initial data (ok.;, oK^) for the system H4-16) and \3.2 1) on a p = po level 


surface are chosen such that they are all smooth on their respective domains, and 
that the system comprised by . i 6| ) and H3.2\) can be guaranteed to be symmetrizable 
hyperbolic in a sufficiently small neighborhood of S^ po in E, as indicated above. Then, 
there exists a unique smooth (local) solution (k;,K^) to d f.lO) and \3.2\) such that 




PO 


= ok;,K* 


l\J?, 


P 0 


= nK 1 


This solution, along with the freely specifiable data 
\f.23fy , uniquely determines n via [ fill ). The fields h y and K i} which can be built 
up from the solution to IffiTfif), H3.2\) and ft4.17\ ) and from the freely specifiable part of 
the data do satisfy the full constraint system in the domain of dependence 

of y po in E. 


To give an explicit example that have immediate physical relevance and where n 
and K \ satisfy globally the inequality n ■ < 0 on some initial data surface restrict 
our considerations to the space of near Schwarzschild configurations. Recall that the 
Schwarzschild metric in the Kerr-Schild form on a Minkowski background (R 4 ,^ a &) 
read as 

gab = Vab + 27744 , (4.29) 

where H = Mjr , which (apart from the singularity at the center, r = 0) is smooth 
and 4 is a null vector field with respect to both g ab and rj ab . Choose now a t = const 
Kerr-Schild time slice as our initial data surface E in a generic Kerr-Schild setup, 
and assumed that the freely specifiable part of the data approximates that of the 
Schwarzschild solution and E is foliated by “r = const ” level surfaces. In particular, 
if in addition the relation — « 0 is satisfied on E, then, as it was shown in [16], for 
these “near Schwarzschild configurations” the relation 


K \ 2(1 + 277) 


K 


1 + 77 

holds on any of the t = const Kerr-Schild time slices. 


(4.30) 
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5 Final remarks 


In this paper the constraints for smooth [n+ 1]-dimensional (with n > 3) Riemannian 
or Lorentzian spaces satisfying the Einstein equations were considered. The initial 
data surface was assumed to be foliated by a one-parameter family of codimension- 
one surfaces, which, for simplicity, were assumed to be orientable either compact 
without boundary in £ or non-compact and infinite. It was shown, regardless whether 
the primary space was Riemannian or Lorentzian, that the constraints can be put 
into the form of a coupled hyperbolic-parabolic system (if the mean curvature of 
the foliating surfaces is positive) or, alternatively, into a Friedrichs symmetrizable 
hyperbolic system and a subsidiary algebraic relation. In both cases those conditions 
which guarantee the (local) existence and uniqueness of smooth solutions to these 
evolutionary systems were identified. 

Notably both of the above proposed evolutionary methods provide immediate 
answer to the first part of the questions raised in | 20 j (see also the second paragraph 
of the introduction of this paper) by providing a distinguished splitting of the basic 
variables h t j and K,j , or that of IV, iV\ 7 ^, k, k,:, K/, Kjj. For instance, in case of 
the hyperbolic-parabolic system N, k, and are constrained whereas N l , 7 y, k, and 

O ^ . 

Kij are freely specifiable throughout £. The initial data (oiV, okj, 0K7), with qN > 
0 , is also freely specifiable on one of the level surfaces of the foliation S? p of £. 
Alternatively, in case of the symmetrizable hyperbolic system the variables k. k, and 
K l [ are constrained whereas N, N l , 7 ^ and K^ are freely specifiable throughout £. 
Note that the initial data for (14. 16j) and (13.21) is also at our disposal, i.e. the values 
of ok.; and oK \ are freely specifiable on one of the level surfaces of the foliation 5? p 
of £. 

Note that the method built on the first order symmetrizable hyperbolic system 
(I4.16P and (13. 2 1) appears to be applicable in deducing initial data for black hole con¬ 
figurations PS]. One of the advantages of using the pertinent evolutionary setup is 
that—as opposed to the requirements of the conformal method where one needs to 
specify boundary data in the asymptotic region, as well as, in the strong field regime 
inside the black holes—it suffices to specify initial data only on a distant sphere sur¬ 
rounding the binary system in the asymptotic region. It is expected that the use of the 
proposed evolutionary approach will yield a reduction of the spurious gravitational 
wave content of numerical simulations. 

Note that as the momentum constraint takes the form of a first order symmetric 
hyperbolic system—see also analogous systems in mm — the following interesting 
conceptual issue is raised. Namely, hyperbolicity in the conventional setup is always 
related to the causal structure determined by a, possibly evolving, Lorentzian space- 
time metric. For the first glance the hyperbolicity of the momentum constraint, as 
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given by (13. 1 j) - (I3.2D . appears to be significantly different as—regardless of the signa¬ 
ture of the primary space—the metric induced on £ is simply Riemannian that has 
nothing to do with causality. 

Notice, however, that by assuming the existence of a foliation S? p of £ we also 
introduced a normal vector field h l that is well-defined (apart from centers) on £. 
Taking then into account the well-known result that with the help of such a vector 
field it is always possible to define a Lorentzian metric h %J on £ via the relation (see, 
e.g. PH) 

hij = hij - (1 + a) fiifij , (5.1) 

where a is a positive real function on £. (In particular, if a = 1 the vector field n l 
gets to be a unit norm timelike vector field with respect to the Lorentzian metric hij.) 

Note that in virtue of (12.11) and (13.61) . the latter determining the characteristic cone 
for the system fl3.1j) - fl3.2D , a Lorentzian metric h r j , with a = n — 1 > 2 (for n > 3), 
entered implicitly into our discussions. This provides then an immediate explanation 
of the appearance of hyperbolicity in context of the momentum constraint. 

The linearity of the Hamiltonian constraint (13.81) in n may play important role 
in canonical quantization of gravity. For instance, in Dirac’s approach that is known 
to apply to constrained systems such as Einstein’s theory of gravity, the constraints 
are ignored at the classical level and the physical quantum states are supposed to be 
determined by solving the quantum constraints on a kinematical Hilbert space. It is 
of obvious interest to know whether the linearity of the Hamiltonian constraint in n 
could simplify the identification of the true physical quantum states. 

As already indicated in section EJ it would be of obvious interest to identify those 
conditions which could guarantee the global existence and uniqueness—and, possibly, 
the asymptotically Euclidean character and/or the regularity at centers—of solutions 
to the evolutionary systems introduced in this paper. 

The clear up of all the above mentioned open problems, in particular, those in the 
last two paragraphs—and the involvement of foliations the level surfaces of which are 
compact with boundary in £ (see also footnote El)— would certainly deserve further 
attention. 
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